Bandgap Calculation and Experimental Analysis of Piezoelectric Phononic Crystals Based on Partial Differential Equations

Focusing on the bending wave characteristic of plate–shell structures, this paper derives the complex band curve of piezoelectric phononic crystal based on the equilibrium differential equation in the plane stress state using COMSOL PDE 6.2. To ascertain the computational model’s accuracy, the computed complex band curve is then cross-validated against real band curves obtained through coupling simulations. Utilizing this model, this paper investigates the impact of structural and electrical parameters on the bandgap range and the attenuation coefficient in the bandgap. Results indicate that the larger surface areas of the piezoelectric sheet correspond to lower center bands in the bandgap, while increased thickness widens the attenuation coefficient range with increased peak values. Furthermore, the influence of inductance on the bandgap conforms to the variation law of the electrical LC resonance frequency, and increased resistance widens the attenuation coefficient range albeit with decreased peak values. The incorporation of negative capacitance significantly expands the low-frequency bandgap range. Visualized through vibration transfer simulations, the vibration-damping ability of the piezoelectric phononic crystal is demonstrated. Experimentally, this paper finds that two propagation modes of bending waves (symmetric and anti-symmetric) result in variable voltage amplitudes, and the average vibration of the system decreases by 4–5 dB within the range of 1710–1990 Hz. The comparison between experimental and model-generated data confirms the accuracy of the attenuation coefficient calculation model. This convergence between experimental and computational results emphasizes the validity and usefulness of the proposed model, and this paper provides theoretical support for the application of piezoelectric phononic crystals in the field of plate–shell vibration reduction.


Introduction
Vibration is an inherent phenomenon that occurs during the operation of mechanical systems.While beneficial vibrations can be harnessed for various purposes, detrimental vibrations may cause structural damage, disrupt system functionality, and significantly shorten the lifespan of mechanical systems.Thus, effective vibration control is paramount in ensuring the stability, safety, and operational comfort of mechanical structures.
Vibration isolation can be mainly divided into active vibration isolation and passive vibration isolation [1][2][3].Active vibration isolation involves the use of closed-loop control systems, utilizing actuators and other power output components to counteract vibrations.In contrast, passive vibration isolation primarily relies on structural arrangements between the vibration source and the isolated object to attenuate vibration transmission.Phononic crystals, as an innovative type of passive vibration damping material, consist of periodic structures made from two or more media.These structures modulate elastic waves at specific frequencies, creating a frequency range known as the forbidden band or bandgap [4,5].By employing a thoughtful structural design, phononic crystals can effectively damp vibrations within designated frequency bands, introducing a novel approach to passive vibration isolation.The bandgap range of traditional phononic crystals primarily depends on the unit size and material composition.Once the material and structure are set in actual production, the elastic wave bandgap remains fixed.Furthermore, traditional phononic crystals require substantial additional mass, making it difficult to modify the original structure or to use them in low-strength structures, limiting their practical applicability.Consequently, novel metamaterials like piezoelectric phononic crystals are attracting increasing attention from researchers.
Piezoelectric phononic crystals utilize piezoelectric materials as scatterers.These smart materials can convert mechanical energy into electrical energy and dissipate it through shunt damping in external circuits [6,7].This process modulates elastic waves via electrical resonance and active control without altering the original structure [8,9].Li Zhengyang [10] systematically summarized the current research status of smart piezoelectric materials and classified them into single-type piezoelectric phononic crystals and metamaterials, embedded piezoelectric/elastic composite-type piezoelectric phononic crystals and metamaterials, and composite piezoelectric phononic crystals and metamaterials consisting of externally attached piezoelectric sheets and elastic structures.Patch-type piezoelectric sheets have a small additional mass, facilitating the implementation of external circuits and active control methods.Yi [11] categorized piezoelectric phononic crystals into shuntdamped and sensing-actuated types based on their working principles.The shunt-damped type utilizes the intrinsic capacitance of the piezoelectric sheet to create an electrical resonance effect with external circuits, thereby maximizing the attenuation of elastic waves.Wang [12] derived the transfer function of the piezoelectric phononic crystal semi-active control strategy using a digital control circuit and proposed a "pole-zero" control strategy; Sheng-Bing [13] used an anti-symmetric arrangement of piezoelectric sheets to control the wave propagation and attenuation characteristics of flexible beams and explored the changes in their vibration isolation performance.
Researchers globally have extensively studied piezoelectric phononic crystals.Shuntdamped piezoelectric phononic crystals use shunt circuit damping as the primary vibration isolation mechanism with the bandgap range and attenuation capability within the bandgap serving as key metrics for assessing their vibration isolation performance.Currently, commonly employed bandgap calculation methods include the Transfer Matrix Method [14], Spectral Element Method, Plane Wave Expansion Method [15], and Finite Element Method [16][17][18][19].Kargozarfard [20] proposed an enhanced plane-wave expansion method that addresses the limitations of conventional plane-wave techniques and enables precise control of the bandgap structure via spatiotemporal modulation characteristics.Shen [21] utilized the "negative equivalent elastic modulus" property of piezoelectric sheets to approximate the local resonance bandgap range and verified it numerically.Panxue [22] derived the vibration transmission rate of electromechanical coupled systems under finite periods based on Timoshenko Beam Theory and investigated the effects of structural and electrical parameters on vibration damping performance.Zhiwen [23] developed a dispersion calculation model for piezoelectric phononic crystal plates using Mindlin's theory and the plane wave expansion method, calculating shear correction coefficients via transcendental equations.Sugino [24] formulated the governing equations of a twodimensional fully coupled electromechanical system and deduced the limiting-state form of the edge frequencies of the bandgap.Shengbing [25,26] created an acoustic transmission loss coefficient model based on dynamic equivalence theory.Milić [27] introduced a finite element analysis method founded on NURBS and Reissner-Mindlin kinematics to address the behavior of active composite laminates exhibiting electro-mechanical coupling effects.Nestorović [28] proposed an ABAQUS-Python co-simulation method that efficiently handles complex thin-walled structures and mitigates the locking effect.Since piezoelectric phononic crystals encompass both mechanical and electrical fields, coupling multi-physical fields is essential for simulation calculations.Analyzing the attenuation coefficient within the bandgap is crucial due to the damping properties of the external shunt circuit.Thus, it is necessary to analyze the piezoelectric phononic crystal using complex band theory.Current bandgap simulations for piezoelectric phononic crystals primarily rely on the Kirchhoff plate assumption to derive finite element models, which involve extensive mathematical operations and present challenges with solving eigenvalue functions.
In summary, this paper derives the formula for the complex band calculation of piezoelectric phononic crystals based on the equilibrium differential equation under plane stress conditions and the intrinsic equation of the piezoelectric sheet.By introducing the Bloch boundary condition through parameter changes, it utilizes the strong coupling capability of the mathematical partial differential equation (PDE) module to quickly calculate the bandgap range and attenuation coefficient magnitude within the bandgap.The obtained bandgap range is verified through the multi-physics coupling of COMSOL's structural mechanics and electrical modules.The effects of different inductance, resistance, and capacitance parameters on the bandgap range and attenuation coefficient in the bandgap are analyzed for various coverage and external shunt circuits, and the resulting influence laws are summarized.The impact of electrical parameters on vibration isolation performance is validated by comparing the simulation of vibration transfer curves from finite element software with experimental results.
The research demonstrates that the complex energy band calculation method utilized in this paper exhibits excellent computability for determining the bandgap position, range, and attenuation coefficient in piezoelectric phonon crystals.Furthermore, the modeling with COMSOL finite element simulation software enhances the calculation method's applicability to piezoelectric phonon crystals on complex patch surfaces, such as acoustic black hole surfaces.It is proven that piezoelectric phononic crystals exhibit excellent vibration damping performance under finite period arrangements.Through experiments, we analyze the signals generated by piezoelectric plates in different connection modes and investigate their underlying causes.The accuracy of the attenuation coefficient is validated by comparing it with the vibration transmission curve measured in experiments.

Calculated Parameters for the Complex Band of Piezoelectric Phononic Crystals
As illustrated in Figure 1, the 2D piezoelectric phononic crystal plate primarily consists of a substrate plate, a piezoelectric ceramic sheet (PZT), and an external shunt circuit.Figure 2 shows the schematic diagram of its unit cell, where a thin plate functions as the substrate, and the piezoelectric sheet along with the external shunt circuit constitutes the scatterer.The piezoelectric ceramic sheet is symmetrically attached to the upper and lower surfaces of the substrate plate and is simultaneously connected to the external circuit.The width and thickness of the substrate plate for a single unit cell are denoted by l b and h b , respectively, while those of the piezoelectric sheet are denoted l p and h p .The complex impedance of the external shunt circuit is denoted by Z p .
In the force-electric coupling system of piezoelectric phononic crystals, the piezoelectric sheet, as a smart material, can convert mechanical energy into electrical energy.When the substrate plate vibrates at an angular frequency ω with simple harmonic waves, the piezoelectric sheet on its upper and lower surfaces uses the positive piezoelectric effect to convert mechanical energy into electrical energy.This energy is dissipated through the external shunt damping circuit.According to the structural properties of the phononic crystal, which possesses a two-dimensional periodic structure, the displacement field can be decomposed by applying Bloch's theorem as follows:  In the force-electric coupling system of piezoelectric phononic crystals, the piezoelectric sheet, as a smart material, can convert mechanical energy into electrical energy.When the substrate plate vibrates at an angular frequency ω with simple harmonic waves, the piezoelectric sheet on its upper and lower surfaces uses the positive piezoelectric effect to convert mechanical energy into electrical energy.This energy is dissipated through the external shunt damping circuit.According to the structural properties of the phononic crystal, which possesses a two-dimensional periodic structure, the displacement field can be decomposed by applying Bloch's theorem as follows: For analyzing real bands in COMSOL's structural mechanics mode, the coupling between multiphysics fields is utilized by setting up structural mechanics, electrostatics, and circuit modules.The piezoelectric effect facilitates state coupling between the mechanical and electrical fields.By scanning the wavevector k, different eigenmodes of the structure are obtained.When no eigenmode exists in a corresponding interval, this interval is identified as the bandgap.This method is known as real band analysis or ω(k)-analysis.In this method, only real numbers are used for scanning the wavevector k.This allows for determining the position of the bandgap in the energy band analysis, but it does not provide the specific attenuation coefficient within the bandgap.For piezoelectric phononic crystals under force-electrical coupling, real band analysis cannot reflect the width of the local resonance bandgap or the damping effect in the shunt circuit, only the center of the local resonance bandgap.Consequently, real band analysis cannot determine the range of the local resonance bandgap or the attenuation capability of elastic waves within this range.
Obtaining the eigenwave vector k by specifying the angular frequency ω is known as complex band analysis or k(ω) analysis.This method yields a complex expression of the eigenwave vector:  In the force-electric coupling system of piezoelectric phononic crystals, the piezoelectric sheet, as a smart material, can convert mechanical energy into electrical energy.When the substrate plate vibrates at an angular frequency ω with simple harmonic waves, the piezoelectric sheet on its upper and lower surfaces uses the positive piezoelectric effect to convert mechanical energy into electrical energy.This energy is dissipated through the external shunt damping circuit.According to the structural properties of the phononic crystal, which possesses a two-dimensional periodic structure, the displacement field can be decomposed by applying Bloch's theorem as follows: For analyzing real bands in COMSOL's structural mechanics mode, the coupling between multiphysics fields is utilized by setting up structural mechanics, electrostatics, and circuit modules.The piezoelectric effect facilitates state coupling between the mechanical and electrical fields.By scanning the wavevector k, different eigenmodes of the structure are obtained.When no eigenmode exists in a corresponding interval, this interval is identified as the bandgap.This method is known as real band analysis or ω(k)-analysis.In this method, only real numbers are used for scanning the wavevector k.This allows for determining the position of the bandgap in the energy band analysis, but it does not provide the specific attenuation coefficient within the bandgap.For piezoelectric phononic crystals under force-electrical coupling, real band analysis cannot reflect the width of the local resonance bandgap or the damping effect in the shunt circuit, only the center of the local resonance bandgap.Consequently, real band analysis cannot determine the range of the local resonance bandgap or the attenuation capability of elastic waves within this range.
Obtaining the eigenwave vector k by specifying the angular frequency ω is known as complex band analysis or k(ω) analysis.This method yields a complex expression of the eigenwave vector: For analyzing real bands in COMSOL's structural mechanics mode, the coupling between multiphysics fields is utilized by setting up structural mechanics, electrostatics, and circuit modules.The piezoelectric effect facilitates state coupling between the mechanical and electrical fields.By scanning the wavevector k, different eigenmodes of the structure are obtained.When no eigenmode exists in a corresponding interval, this interval is identified as the bandgap.This method is known as real band analysis or ω(k)-analysis.In this method, only real numbers are used for scanning the wavevector k.This allows for determining the position of the bandgap in the energy band analysis, but it does not provide the specific attenuation coefficient within the bandgap.For piezoelectric phononic crystals under force-electrical coupling, real band analysis cannot reflect the width of the local resonance bandgap or the damping effect in the shunt circuit, only the center of the local resonance bandgap.Consequently, real band analysis cannot determine the range of the local resonance bandgap or the attenuation capability of elastic waves within this range.
Obtaining the eigenwave vector k by specifying the angular frequency ω is known as complex band analysis or k(ω) analysis.This method yields a complex expression of the eigenwave vector: where k r and k i represent the real and imaginary parts of the wave vector, respectively.Compared to real band analysis, this method introduces the imaginary part of the wave vector to describe the amplitude changes in the vibration of the elastic wave.For a forceelectric coupled system analyzed using complex bands, the intrinsic wave vector k obtained by scanning the angular frequency can reflect the range of the local resonance bandgap due to circuit resonance as well as the magnitude of the attenuation coefficient within the bandgap.
For the piezoelectric phononic crystal analyzed in this paper, the thicknesses of the substrate plate and the piezoelectric sheet are much smaller than the minimum size of the midplane.Therefore, the analysis can be simplified using thin-plate theory: where Y is the Young's modulus of the substrate plate, µ is the Poisson's ratio of the substrate plate, u 1 is the displacement component in the x-direction, and u 2 is the displacement component in the y-direction.The equilibrium differential equation of any point in the plane stress state is Substituting Equation ( 3) into (4) yields By introducing a parametric transformation equation, the Bloch boundary condition is converted into a general continuum periodic boundary condition.This transformation eliminates the effect of the scanning dependent variable, specifically the intrinsic wave vector, in the boundary condition.The displacement field parametric transformation equation can be expressed as [29] u(r, t) = u k (r)e ik(cos θx+sin θy)−iωt (6) The parametrically transformed displacement field is substituted into Equation ( 5) and simplified, as detailed in Appendix A. The piezoelectric sheets are polarized along the z-axis with the assumption that they are uniformly polarized.The dimension in the polarization direction, i.e., the thickness of the piezoelectric sheet, is significantly smaller than its length and width.The first type of piezoelectric equation system can be expressed as where s E 11 , s E 12 , and s E 66 represent the flexibility coefficient of the piezoelectric sheet under constant electric field strength, d 31 represents the piezoelectric constant, and ε T 33 represents the dielectric constant under constant stress.
Assuming the external excitation wavelength is much larger than the length and width of the piezoelectric sheet, placing the piezoelectric sheet within the approximate wavelength range, the electric field strength E 3 can be considered uniform throughout.As shown in Figure 2, the external circuit connection can be modeled as an AC power supply and capacitance.The upper and lower surfaces of the piezoelectric sheet are connected in parallel, and the voltage across the two terminals of the piezoelectric sheet is equal to the superimposed charge.Based on the external circuit connection and the voltage-current relationship, the following can be derived: and The combinations yield Substituting the equation for the potential shift in terms of electric field strength into the first type of piezoelectric system equations yields the following: where C T p = ε T 33 A p /h p is the intrinsic capacitance of the piezoelectric sheet under mechanically free boundary conditions, corresponding to the constant stress state.To analyze the local resonance bandgap generation mechanism of piezoelectric phononic crystals under a force-electric coupling system, it is essential to consider the intrinsic capacitance formed by the piezoelectric sheet and the inductance of the external circuit, which together form an electrical resonance under an AC source.Therefore, studying the intrinsic capacitance of the piezoelectric sheet is crucial for determining the center band of the local resonance bandgap theoretically.The capacitance of the piezoelectric sheet is dynamic during operation, and it is always situated between the two limiting states: the completely free and completely cut-off states.Consequently, its dynamic capacitance lies between the capacitance of the constant stress state and the constant strain state.The boundary conditions determined by the first type of piezoelectric equations represent the completely free boundary in the constant stress state.Similarly, the complete cut-off state in the constant strain state is calculated by the third type of piezoelectric equations.The piezoelectric sheet capacitance in these states is given by The dynamic capacitance of the piezoelectric sheet always lies within the two limit state intervals.The presence of dynamic capacitance during operation can be corrected for the local resonance bandgap center position by introducing a capacitance correction factor τ. By substituting the electric field strength into the stress-strain equations for coupled mechanical and electrical field modeling, the first type of piezoelectric equation set can be transformed into where Similarly, the control differential equations in the piezoelectric domain can be obtained by substituting the parameter-transformed displacement field expressions.These equations can then be combined with the control differential equations in the substrate plate domain.For the force-electrically coupled piezoelectric phononic crystal, the COMSOL PDE module can be used to quickly solve the second-order partial differential equations and impose general continuity conditions.The overall control equations can be listed as follows: where For substrate panels: For piezoelectric sheets:

Model Validation
In this paper, the COMSOL solid mechanics, electrostatics, and circuits modules are used to perform real band calculations for piezoelectric phononic crystals through the piezoelectric effect coupled with multiphysics fields.The substrate plate is made of epoxy resin, and the piezoelectric sheet model is PZT-5H.The specific material parameters are shown in Tables 1 and 2. The external shunt circuit employs an LR shunt damping circuit, integrating structural mechanics with the electrical module.Periodic boundary conditions are implemented through structural mechanics, as illustrated in Figure 3.The circuit terminal is positioned on the upper and lower surfaces of the piezoelectric plate via the electrostatic module, and it is coupled with the electrical device through the circuit module to resolve the real band.The solution process includes some in-plane modes, but the generation of positiv and negative piezoelectric effects and the realization of vibration isolation occur in th out-of-plane modes.Since the piezoelectric phononic crystal provides vibration isolatio for bending waves (out-of-plane modes), it is necessary to establish a screening criterio for bending waves and exclude the in-plane modes [30]:   (17 where, u, v, and w, respectively, represent the displacement of the body element at an point in the x, y, and z direction; u * , v * , and w * , respectively, represent the body element a any point in the direction of the displacement of the conjugate complex of the x, y, z direc tions, taking β_Modal greater than or equal to 0.9 as the out-of-plane mode.
As shown in Figure 4a, due to the spatial symmetry of the lattice structure of phono crystals, for two-dimensional plane problems, the band curve can be obtained by simpl scanning the wave vector k along the boundary of the minimum irreducible Brilloui zone, as illustrated in Figure 4b.Two bandgap ranges are generated in the X Γ direction The first bandgap (LG1, 581-590 Hz) is dominated by the local resonance effect, resultin from the electrical resonance between the intrinsic capacitance of the piezoelectric shee and the inductance of the external circuit.This bandgap is extremely narrow, because rea band analysis cannot reflect the resistive damping effects of the piezoelectric phononi crystal.The second bandgap (BG1, 730-990 Hz) is dominated by the Bragg scattering e fect, where the elastic wave is confined by the periodic medium, resulting in attenuation The MΓ direction generates a bandgap identical in mechanism and position to the firs bandgap in the X Γ direction.This is mainly because the local resonance bandga formed by electrical resonance is not directionally selective and can form an omnidirec tional bandgap.The solution process includes some in-plane modes, but the generation of positive and negative piezoelectric effects and the realization of vibration isolation occur in the out-of-plane modes.Since the piezoelectric phononic crystal provides vibration isolation for bending waves (out-of-plane modes), it is necessary to establish a screening criterion for bending waves and exclude the in-plane modes [30]: where, u, v, and w, respectively, represent the displacement of the body element at any point in the x, y, and z direction; u * , v * , and w * , respectively, represent the body element at any point in the direction of the displacement of the conjugate complex of the x, y, z directions, taking β_Modal greater than or equal to 0.9 as the out-of-plane mode.
As shown in Figure 4a, due to the spatial symmetry of the lattice structure of phonon crystals, for two-dimensional plane problems, the band curve can be obtained by simply scanning the wave vector k along the boundary of the minimum irreducible Brillouin zone, as illustrated in Figure 4b.Two bandgap ranges are generated in the ΓX direction.The first bandgap (LG1, 581-590 Hz) is dominated by the local resonance effect, resulting from the electrical resonance between the intrinsic capacitance of the piezoelectric sheet and the inductance of the external circuit.This bandgap is extremely narrow, because real band analysis cannot reflect the resistive damping effects of the piezoelectric phononic crystal.The second bandgap (BG1, 730-990 Hz) is dominated by the Bragg scattering effect, where the elastic wave is confined by the periodic medium, resulting in attenuation.The MΓ direction generates a bandgap identical in mechanism and position to the first bandgap in the ΓX direction.This is mainly because the local resonance bandgap formed by electrical resonance is not directionally selective and can form an omnidirectional bandgap.
Because the real band cannot reflect the impact of resistive damping, the resistance R is set to 0 Ω during the complex band analysis to maintain consistency between the real and complex bands.For the complex band analysis, the elastodynamic fluctuation equations and the electric field are coupled through the COMSOL PDE module using second-order partial differential equations.The wave equation is input in the form of Equation ( 16), and the boundary is constrained by Dirichlet boundary conditions (general continuity boundaries).Because the real band cannot reflect the impact of resistive damping, the resistance R is set to 0Ω during the complex band analysis to maintain consistency between the real and complex bands.For the complex band analysis, the elastodynamic fluctuation equations and the electric field are coupled through the COMSOL PDE module using secondorder partial differential equations.The wave equation is input in the form of Equation ( 16), and the boundary is constrained by Dirichlet boundary conditions (general continuity boundaries). -- Similarly, the in-plane modes are screened out by data screening, and the eigenwave vectors of the out-of-plane modes are solved for k.When the eigenwave vectors are all real, the wave can propagate without attenuation.When the eigenwave vectors contain an imaginary part, it indicates that the wave will attenuate during transmission.The propagation coefficient is defined as where a is the lattice constant.The real part represents the phase change in the wave vector, and the imaginary part indicates the magnitude of the wave attenuation coefficient during propagation.The phase coefficient changes abruptly within the local resonance bandgap, while the phase constant remains unchanged within the Bragg scattering bandgap.By comparing the real band with the center of the local resonance bandgap of the complex band, the dynamic capacitance correction coefficient is determined to be 0.85.The corrected complex band curve is plotted as shown in Figure 5. Similarly, the in-plane modes are screened out by data screening, and the eigenwave vectors of the out-of-plane modes are solved for k.When the eigenwave vectors are all real, the wave can propagate without attenuation.When the eigenwave vectors contain an imaginary part, it indicates that the wave will attenuate during transmission.The propagation coefficient is defined as where a is the lattice constant.The real part represents the phase change in the wave vector, and the imaginary part indicates the magnitude of the wave attenuation coefficient during propagation.The phase coefficient changes abruptly within the local resonance bandgap, while the phase constant remains unchanged within the Bragg scattering bandgap.By comparing the real band with the center of the local resonance bandgap of the complex band, the dynamic capacitance correction coefficient is determined to be 0.85.The corrected complex band curve is plotted as shown in Figure 5.As shown in Figure 5, the phase coefficients in the 0 • complex band curves calculated by the PDE module exhibit a high degree of similarity with the real band curves obtained from ΓX directional scanning using the structural mechanics coupling module.This consistency verifies the accuracy of the complex band calculation method used in this paper.Since the local resonance bandgap generated by piezoelectric sheet force-electric coupling is not directionally selective, and this paper focuses on the local resonance bandgap generated by force-electric coupling, the simulation of the PDE module in this paper is based on the local resonance bandgap at 0 • to simplify the model calculation.The position of the local resonance bandgap and the decay coefficient inside the bandgap are analyzed under the condition of 0 • directional scanning.As shown in Figure 5, the phase coefficients in the 0° complex band curves calculated by the PDE module exhibit a high degree of similarity with the real band curves obtained from ΓΧ directional scanning using the structural mechanics coupling module.This consistency verifies the accuracy of the complex band calculation method used in this pa per.Since the local resonance bandgap generated by piezoelectric sheet force-electric cou pling is not directionally selective, and this paper focuses on the local resonance bandgap generated by force-electric coupling, the simulation of the PDE module in this paper is based on the local resonance bandgap at 0° to simplify the model calculation.The position of the local resonance bandgap and the decay coefficient inside the bandgap are analyzed under the condition of 0° directional scanning.

Effect of Coverage on Bandgap
According to the bandgap generation principle of piezoelectric phononic crystals and the series-parallel circuit resonance equation, the coverage is directly related to the posi tion of the local resonance bandgap.The coverage is defined as In order to exclude the influence produced by the Bragg bandgap, this paper selects the substrate plate material as the aluminum plate material close to the elastic modulus and Poisson's ratio of PZT-5H, the lattice constant (width of the substrate plate, lb) is 80 mm, the thickness of the substrate plate, hb, is 4 mm, and the intrinsic capacitance of the piezoelectric sheet can be expressed in terms of the coverage ratio as

Effect of Coverage on Bandgap
According to the bandgap generation principle of piezoelectric phononic crystals and the series-parallel circuit resonance equation, the coverage is directly related to the position of the local resonance bandgap.The coverage is defined as In order to exclude the influence produced by the Bragg bandgap, this paper selects the substrate plate material as the aluminum plate material close to the elastic modulus and Poisson's ratio of PZT-5H, the lattice constant (width of the substrate plate, l b ) is 80 mm, the thickness of the substrate plate, h b , is 4 mm, and the intrinsic capacitance of the piezoelectric sheet can be expressed in terms of the coverage ratio as The resonance center interval, which is the location of the center band of the local resonance bandgap, can be obtained as When the thickness of the piezoelectric sheet is fixed and used with the external shunt circuit, the resonance interval becomes a power exponential function of the coverage.The thickness of the piezoelectric sheet is set to 1 mm, the inductance of the external circuit is 0.6 H, the resistance is 30 Ω, and the coverage ranges from 0.4 to 0.9.The complex band attenuation coefficient is calculated at intervals of 0.1.The resulting attenuation coefficient curve is shown in Figure 6.
When the thickness of the piezoelectric sheet is fixed and used with the external shun circuit, the resonance interval becomes a power exponential function of the coverage.Th thickness of the piezoelectric sheet is set to 1 mm, the inductance of the external circuit i 0.6 H, the resistance is 30 Ω, and the coverage ranges from 0.4 to 0.9.The complex band attenuation coefficient is calculated at intervals of 0.1.The resulting attenuation coefficien curve is shown in Figure 6.As shown in Figure 6, changes in coverage primarily affect the position of the cente band of the local resonance bandgap with minimal impact on the attenuation coefficien in the forbidden band interval.The effect on the center band position is mainly due to th change in the intrinsic capacitance of the piezoelectric sheet caused by the variation in th coverage ratio, which leads to changes in the resonance interval.To further explore th effect of coverage and eliminate the influence of changes in the intrinsic capacitance of th piezoelectric sheet on the results, normalized coverage is used: where * p h is the thickness of the piezoelectric sheet after the change, and hp is 1 mm.A this time, the position of the center band of the local resonance bandgap is As shown in Figure 6, changes in coverage primarily affect the position of the center band of the local resonance bandgap with minimal impact on the attenuation coefficient in the forbidden band interval.The effect on the center band position is mainly due to the change in the intrinsic capacitance of the piezoelectric sheet caused by the variation in the coverage ratio, which leads to changes in the resonance interval.To further explore the effect of coverage and eliminate the influence of changes in the intrinsic capacitance of the piezoelectric sheet on the results, normalized coverage is used: where h * p is the thickness of the piezoelectric sheet after the change, and h p is 1 mm.At this time, the position of the center band of the local resonance bandgap is From Equation (24), after applying the normalized coverage setting, the position of the center band of the local resonance bandgap becomes a constant value, excluding the effect of the variation in the intrinsic capacitance of the piezoelectric sheet.
As shown in Figure 7, the position of the center band of the local resonance bandgap is independent of the variation in normalized coverage.As derived from the theory, increasing the normalized coverage results in a wider bandgap and a higher attenuation coefficient within the bandgap.
When the area of the piezoelectric sheet changes but the thickness remains constant, the intrinsic capacitance of the piezoelectric sheet changes, causing the position of the center band of the local resonance bandgap to shift, exhibiting a power-exponential relationship with the coverage rate, while the attenuation coefficients within the forbidden band remain the same.When both the area and thickness of the piezoelectric sheet change, i.e., when the normalized coverage rate changes, the position of the center band of the local resonance bandgap remains unchanged, and the range of the bandgap and the attenuation coefficient within the bandgap increase with the increasing normalized coverage rate.
Materials 2024, 17, x FOR PEER REVIEW 13 of 2 From Equation ( 24), after applying the normalized coverage setting, the position o the center band of the local resonance bandgap becomes a constant value, excluding th effect of the variation in the intrinsic capacitance of the piezoelectric sheet.
As shown in Figure 7, the position of the center band of the local resonance bandga is independent of the variation in normalized coverage.As derived from the theory, in creasing the normalized coverage results in a wider bandgap and a higher attenuatio coefficient within the bandgap.When the area of the piezoelectric sheet changes but the thickness remains constan the intrinsic capacitance of the piezoelectric sheet changes, causing the position of the cen ter band of the local resonance bandgap to shift, exhibiting a power-exponential relation ship with the coverage rate, while the attenuation coefficients within the forbidden ban remain the same.When both the area and thickness of the piezoelectric sheet change, i.e when the normalized coverage rate changes, the position of the center band of the loca resonance bandgap remains unchanged, and the range of the bandgap and the attenuatio coefficient within the bandgap increase with the increasing normalized coverage rate.

Effect of Inductance on Bandgap
According to Equation ( 24), the inductance is a power exponential function of th central band of the local resonance.To further investigate the effect of inductance size o the bandwidth of the elastic wave band and the attenuation coefficient, we choose th width of the substrate plate lb to be 80 mm, the thickness of the substrate plate hb to be mm, the width of the piezoelectric sheet lp to be 60 mm, the thickness of the piezoelectri sheet hp to be 1 mm, and the external shunt circuit to be the LR circuit with a resistance o 30Ω and an inductance range of 0.08-0.30H.The attenuation coefficients for different in ductances are plotted in Figure 8.

Effect of Inductance on Bandgap
According to Equation ( 24), the inductance is a power exponential function of the central band of the local resonance.To further investigate the effect of inductance size on the bandwidth of the elastic wave band and the attenuation coefficient, we choose the width of the substrate plate l b to be 80 mm, the thickness of the substrate plate h b to be 5 mm, the width of the piezoelectric sheet l p to be 60 mm, the thickness of the piezoelectric sheet h p to be 1 mm, and the external shunt circuit to be the LR circuit with a resistance of 30 Ω and an inductance range of 0.08-0.30H.The attenuation coefficients for different inductances are plotted in Figure 8.
From Figure 8, the central band of the local resonance bandgap is a power exponential function of the inductance size.This is mainly because the change in inductance leads to a change in the electrical resonance frequency, thus shifting the central band of the bandgap, which is consistent with the theoretical analysis.As the inductance value increases, the range of the bandgap gradually narrows, and the attenuation coefficient within the bandgap increases slightly.However, the overall fluctuation is not significant, which is mainly because the quality factor of the series resonance circuit, Q = √ L/C/R, increases with the inductance.The overall quality factor of the circuit increases with the inductance and the current: The current response is accelerated, the passband is narrowed, and the resistive energy dissipation interval is reduced, resulting in a narrower local resonance bandgap.From Figure 8, the central band of the local resonance bandgap is a power exponential function of the inductance size.This is mainly because the change in inductance leads to a change in the electrical resonance frequency, thus shifting the central band of the bandgap, which is consistent with the theoretical analysis.As the inductance value increases, the range of the bandgap gradually narrows, and the attenuation coefficient within the bandgap increases slightly.However, the overall fluctuation is not significant, which is mainly because the quality factor of the series resonance circuit, increases with the inductance.The overall quality factor of the circuit increases with the inductance and the current:

Effect of Resistance on Bandgap
From the previous section, it is clear that changes in resistance do not affect the location of the center band of the local resonance bandgap.However, the main vibration isolation unit of the system is damped vibration isolation, and the energy dissipation method is mainly related to the resistance.The inductance of the external circuit is set to 0.2 H, and the resistance is varied from 0 to 300 ohms, with the size parameters being the same as in the previous section.The system attenuation coefficient is plotted accordingly.
From Figure 9, the larger the resistance, the wider the local resonance bandgap range, but the attenuation coefficient within the bandgap decreases.This is mainly because the circuit forms an electrical resonance with the intrinsic capacitance through the external circuit inductance and dissipates vibrational energy through the resistance.An increase in resistance results in a decrease in the peak current at resonance and a decrease in the peak attenuation coefficient within the bandgap.However, an increase in resistance can also significantly enhance the energy dissipation of the system, allowing it to dissipate energy at smaller currents, thereby widening the bandgap range.
tion of the center band of the local resonance bandgap.However, the main vibration isolation unit of the system is damped vibration isolation, and the energy dissipation method is mainly related to the resistance.The inductance of the external circuit is set to 0.2 H, and the resistance is varied from 0 to 300 ohms, with the size parameters being the same as in the previous section.The system attenuation coefficient is plotted accordingly.
From Figure 9, the larger the resistance, the wider the local resonance bandgap range, but the attenuation coefficient within the bandgap decreases.This is mainly because the circuit forms an electrical resonance with the intrinsic capacitance through the external circuit inductance and dissipates vibrational energy through the resistance.An increase in resistance results in a decrease in the peak current at resonance and a decrease in the peak attenuation coefficient within the bandgap.However, an increase in resistance can also significantly enhance the energy dissipation of the system, allowing it to dissipate energy at smaller currents, thereby widening the bandgap range.

Effect of Generalized Capacitance on Bandgap
The electromechanical coupling coefficient is the geometric mean of the reciprocal elastic dielectric energy density over the elastic self-energy density and the dielectric selfenergy density, indicating the efficiency of conversion between mechanical and electrical energy.Limited by the electromechanical coupling coefficient, the width of the local resonance bandgap formed by the LR shunt circuit is narrow.According to the operating con-

Effect of Generalized Capacitance on Bandgap
The electromechanical coupling coefficient is the geometric mean of the reciprocal elastic dielectric energy density over the elastic self-energy density and the dielectric selfenergy density, indicating the efficiency of conversion between mechanical and electrical energy.Limited by the electromechanical coupling coefficient, the width of the local resonance bandgap formed by the LR shunt circuit is narrow.According to the operating conditions and principles in the piezoelectric phononic crystal, the elastic self-energy density of the piezoelectric sheet is The dielectric self-energy density is The mutually elastic dielectric energy density is Its electromechanical coupling coefficient is According to the piezoelectric intrinsic equation, the electromechanical coupling coefficient of the system at this point is obtained: where K a = A p / s E 11 h p represents the stiffness of the system when the piezoelectric elastomer is short-circuited, and C is the intrinsic capacitance of the piezoelectric double wafer.From Equation (30), the electromechanical coupling coefficient is related to the open-circuit stiffness and intrinsic capacitance of the piezoelectric sheet.The value of the electromechanical coupling coefficient generally ranges from 0.3 to 0.7.If the electromechanical coupling coefficient is too large, the system's open-circuit stiffness can change considerably with the electrical boundary conditions, affecting the system's stability.Improving the electromechanical coupling coefficient can be approached from two perspectives: open-circuit stiffness and intrinsic capacitance.When the system is connected in series with negative capacitance, according to the definition of the electromechanical coupling coefficient, it can be determined that where C * represents the system's generalized capacitance.The series negative capacitance does not alter the original open-circuit stiffness and piezoelectric coefficient of the system; however, it modifies the generalized capacitance of the piezoelectric sheet within the system: where α represents the magnitude of the absolute value of the ratio between the negative capacitance and the intrinsic capacitance of the piezoelectric sheet.The introduction of the negative capacitance circuit alters the system's generalized capacitance and thereby shifts the resonance position: The introduction of negative capacitance enhances the electromechanical coupling coefficient and reduces the system's generalized capacitance, thereby effectively broadening the local resonance bandgap.However, this also impacts the position of the electrical resonance, shifting the local resonance bandgap toward higher frequencies.
From Figure 10, the introduction of negative capacitance significantly enhances the bandgap range compared to conventional circuits, owing to the improved electromechanical coupling coefficient.However, this necessitates a larger inductor to compensate for the reduced generalized capacitance, thereby decreasing the overall quality factor of the circuit and leading to the attenuation of the peak attenuation coefficient.

Vibration Transfer Simulation Analysis and Experiment
A five-cycle piezoelectric phononic crystal plate is established, and a finite element model is constructed using COMSOL by coupling the structural mechanics, electrostatics, and circuits modules.The width of the substrate plate in the metacells is 60 mm with a thickness of 3.5 mm, while the piezoelectric sheet measures 40 mm in width and 1 mm in thickness.An LR external shunt circuit is employed, with an inductance of 0.1 H and a resistance of 300 Ω, as shown in Figure 11.

Vibration Transfer Simulation Analysis and Experiment
A five-cycle piezoelectric phononic crystal plate is established, and a finite element model is constructed using COMSOL by coupling the structural mechanics, electrostatics, and circuits modules.The width of the substrate plate in the metacells is 60 mm with a thickness of 3.5 mm, while the piezoelectric sheet measures 40 mm in width and 1 mm in thickness.An LR external shunt circuit is employed, with an inductance of 0.1 H and a resistance of 300 Ω, as shown in Figure 11.
The excitation point is positioned at point A. Acceleration response curves are recorded at point B in the 0 • direction and at point C in the 45 • direction, forming the transfer curve TL as follows: where a out is the output point acceleration and a in is the excitation point acceleration.
Plotting the vibration transfer curve, theoretical calculations place the resonance position at 1736 Hz.As shown in Figure 12, the vibration transfer curve indicates that vibration attenuation occurs at this resonance position in both the 0 • and 45 • directions, thereby validating the theoretical analysis.Among these, the formation diagrams at points A1 and F1 visually illustrate the attenuation process of the bending wave.
A five-cycle piezoelectric phononic crystal plate is established, and a finite element model is constructed using COMSOL by coupling the structural mechanics, electrostatics, and circuits modules.The width of the substrate plate in the metacells is 60 mm with a thickness of 3.5 mm, while the piezoelectric sheet measures 40 mm in width and 1 mm in thickness.An LR external shunt circuit is employed, with an inductance of 0.1 H and a resistance of 300 Ω, as shown in Figure 11.
where aout is the output point acceleration and ain is the excitation point acceleration.Plotting the vibration transfer curve, theoretical calculations place the resonance position at 1736 Hz.As shown in Figure 12, the vibration transfer curve indicates that vibration attenuation occurs at this resonance position in both the 0° and 45° directions, thereby validating the theoretical analysis.Among these, the formation diagrams at points A1 and F1 visually illustrate the attenuation process of the bending wave.The currents through resistors 13 and 15 during the simulation are depicted in Figure 13.Herein, the current peak due to the resonance effect is relatively flat, and resistors No.The through resistors 13 and 15 during the simulation are depicted in Figure 13.Herein, the current peak due to the resonance effect is relatively flat, and resistors No. 13 and No. 15 exhibit the same current peak.
To verify the effectiveness of the piezoelectric phonon crystal in attenuating bending waves, an experimental platform was constructed, as shown in Figure 13.To reduce the complexity of the experiment, the 0 • direction experimental curve, verified in the previous section, is used as the basis for analysis.This experiment employs a B&K data acquisition instrument to output the excitation signal, which is amplified by a power amplifier to drive the electromagnetic actuator.The acceleration signal is collected by a B&K4507B acceleration sensor and processed in the host computer for vibration signal analysis.The experimental test system constructed is depicted in Figure 14.The currents through resistors 13 and 15 during the simulation are depicted in Figure 13.Herein, the current peak due to the resonance effect is relatively flat, and resistors No. 13 and No. 15 exhibit the same current peak.To verify the effectiveness of the piezoelectric phonon crystal in attenuating bending waves, an experimental platform was constructed, as shown in Figure 13.To reduce the complexity of the experiment, the 0° direction experimental curve, verified in the previous section, is used as the basis for analysis.This experiment employs a B&K data acquisition instrument to output the excitation signal, which is amplified by a power amplifier to drive the electromagnetic actuator.The acceleration signal is collected by a B&K4507B acceleration sensor and processed in the host computer for vibration signal analysis.The experimental test system constructed is depicted in Figure 14.This experiment employs a bilateral symmetric paste form.Theoretical analysis suggests that the piezoelectric sheet should produce the same voltage corresponding to the frequency of the excitation signal.However, in practice, the bending wave generated by the excitation process propagates in both symmetric and anti-symmetric modes.The propagation mode of the bending wave affects the amplitude of the voltage generated by different connection forms.The RIGOL DS110Z oscilloscope was employed to measure the voltage values generated by the piezoelectric sheet in various connection forms.The filtered data and their corresponding connection forms are presented below in Figure 15.This experiment employs a bilateral symmetric paste form.Theoretical analysis suggests that the piezoelectric sheet should produce the same voltage corresponding to the frequency of the excitation signal.However, in practice, the bending wave generated by the excitation process propagates in both symmetric and anti-symmetric modes.The propagation mode of the bending wave affects the amplitude of the voltage generated by different connection forms.The RIGOL DS110Z oscilloscope was employed to measure the voltage values generated by the piezoelectric sheet in various connection forms.The filtered data and their corresponding connection forms are presented below in Figure 15.
The frequency of the voltage generated by all three connection forms aligns with the frequency of the excitation source.However, the amplitude of the voltage produced by UCF is approximately the sum of the voltages produced by the PPF and PNF, as voltage can be produced in a unilateral mode regardless of whether the wave propagates in a symmetrical or anti-symmetrical mode.When the voltage amplitude is small, interference from the clutter signal dominates, resulting in a weak resonance phenomenon and drastically reduced vibration isolation effect.Thus, the PPF is not used in this experiment.The intrinsic capacitance of the piezoelectric sheet in the unilateral connection form is half that of the PNF, which is not conducive to reducing the resonant frequency.Therefore, this experiment employs the PNF.The external circuit uses a 0.1 H inductor, which has a larger DC internal resistance of approximately 270 Ω.In the 1000-3000 Hz frequency sweep, the attenuation coefficient and vibration transfer curve are shown in Figure 16 below.
The corresponding attenuation coefficient of the interval is 0.1-0.15,and the attenuation of the vibration amplitude after five cycles should be (e −0.1 ~e−0.15 ) 5 times the original one, which is about 5-7 dB, corresponding to the experimentally measured attenuation of 5 dB and verifying the accuracy of the attenuation coefficient values.This experiment employs a bilateral symmetric paste form.Theoretical analysis suggests that the piezoelectric sheet should produce the same voltage corresponding to the frequency of the excitation signal.However, in practice, the bending wave generated by the excitation process propagates in both symmetric and anti-symmetric modes.The propagation mode of the bending wave affects the amplitude of the voltage generated by different connection forms.The RIGOL DS110Z oscilloscope was employed to measure the voltage values generated by the piezoelectric sheet in various connection forms.The filtered data and their corresponding connection forms are presented below in Figure 15.The frequency of the voltage generated by all three connection forms aligns with the frequency of the excitation source.However, the amplitude of the voltage produced by UCF is approximately the sum of the voltages produced by the PPF and PNF, as voltage can be produced in a unilateral mode regardless of whether the wave propagates in a symmetrical or anti-symmetrical mode.When the voltage amplitude is small, interference from the clutter signal dominates, resulting in a weak resonance phenomenon and drastically reduced vibration isolation effect.Thus, the PPF is not used in this experiment.The intrinsic capacitance of the piezoelectric sheet in the unilateral connection form is half that of the PNF, which is not conducive to reducing the resonant frequency.Therefore, this experiment employs the PNF.The external circuit uses a 0.1 H inductor, which has a larger DC internal resistance of approximately 270 Ω.In the 1000-3000 Hz frequency sweep, the attenuation coefficient and vibration transfer curve are shown in Figure 16 below.The corresponding attenuation coefficient of the interval is 0.1-0.15,and the attenuation of the vibration amplitude after five cycles should be (e −0.1 ~e−0.15 ) 5 times the original one, which is about 5-7 dB, corresponding to the experimentally measured attenuation of 5 dB and verifying the accuracy of the attenuation coefficient values.

Conclusions
In this paper, based on the plane stress state, energy band theory, and Bloch's theorem, the partial differential equations of motion of the system under the influence of electrical parameters are derived using wave vector parameter transformations and numerical analysis.The complex band characteristic curves are solved, and the accuracy of the complex band model is verified by comparing them with the real band curves obtained from the multiphysics coupled field model.
Using the attenuation coefficient as an evaluation index, we analyze the effects of various structural and electrical parameters on the attenuation coefficient and examine the coupling influence of each parameter.By simulating vibration transmission, we can more intuitively observe the bending wave propagation process in various scenarios.
Experiments measured the voltage generated by the piezoelectric sheet under different connection forms and analyzed the reasons for its generation.The accuracy of the attenuation coefficient in the complex band model is verified by comparing it with the experimentally measured range and value of the vibration transfer curve.
In this work, we derive a complex energy band model for piezoelectric phononic crystals based on the equilibrium differential equations in the plane stress state using COMSOL PDE.We validate the model's accuracy by comparing it with simulation and experimental data.This paper provides theoretical support for applying piezoelectric phononic crystals in plate-shell vibration reduction, especially for complex patch surfaces.
However, the complex band calculation method used in this paper requires scanning the wave vector directions of 0°, 45°, and 90° to obtain the complete band curve in the twodimensional plane.Due to the extensive computational requirements, only the all-directional bandgap in the 0° direction is analyzed in this study, simplifying the design for

Conclusions
In this paper, based on the plane stress state, energy band theory, and Bloch's theorem, the partial differential equations of motion of the system under the influence of electrical parameters are derived using wave vector parameter transformations and numerical analysis.The complex band characteristic curves are solved, and the accuracy of the complex band model is verified by comparing them with the real band curves obtained from the multiphysics coupled field model.
Using the attenuation coefficient as an evaluation index, we analyze the effects of various structural and electrical parameters on the attenuation coefficient and examine the coupling influence of each parameter.By simulating vibration transmission, we can more intuitively observe the bending wave propagation process in various scenarios.
Experiments measured the voltage generated by the piezoelectric sheet under different connection forms and analyzed the reasons for its generation.The accuracy of the attenuation coefficient in the complex band model is verified by comparing it with the experimentally measured range and value of the vibration transfer curve.
In this work, we derive a complex energy band model for piezoelectric phononic crystals based on the equilibrium differential equations in the plane stress state using COMSOL PDE.We validate the model's accuracy by comparing it with simulation and experimental data.This paper provides theoretical support for applying piezoelectric phononic crystals in plate-shell vibration reduction, especially for complex patch surfaces.
However, the complex band calculation method used in this paper requires scanning the wave vector directions of 0 • , 45 • , and 90 • to obtain the complete band curve in the two-dimensional plane.Due to the extensive computational requirements, only the alldirectional bandgap in the 0 • direction is analyzed in this study, simplifying the design for future work.This approach enables quicker and more accurate acquisition of the complete energy band curve for piezoelectric phononic crystals in a two-dimensional plane.

Figure 2 .
Figure 2. Schematic diagram of the unit cell.

Figure 2 .
Figure 2. Schematic diagram of the unit cell.

Figure 2 .
Figure 2. Schematic diagram of the unit cell.

Figure 4 .
Figure 4. Real band diagram of piezoelectric phonon crystal.(a) Diagram of the 2D irreducible Brillouin zone.(b) Real band curve (shaded area represents the bandgap range).

Figure 4 .
Figure 4. Real band diagram of piezoelectric phonon crystal.(a) Diagram of the 2D irreducible Brillouin zone.(b) Real band curve (shaded area represents the bandgap range).

Figure 8 .
Figure 8. Attenuation coefficient with different inductances.(a) Comparison of discrete values for different inductances.(b) Comparison of continuous values for different inductances (intensity of the color indicates the magnitude of the attenuation coefficient).

Figure 8 .
Figure 8. Attenuation coefficient with different inductances.(a) Comparison of discrete values for different inductances.(b) Comparison of continuous values for different inductances (intensity of the color indicates the magnitude of the attenuation coefficient).

Figure 9 .
Figure 9. Attenuation coefficient with different resistances.(a) Comparison of discrete values for different resistances.(b) Comparison of continuous values for different resistances (intensity of the color indicates the magnitude of the attenuation coefficient).

Figure 9 .
Figure 9. Attenuation coefficient with different resistances.(a) Comparison of discrete values for different resistances.(b) Comparison of continuous values for different resistances (intensity of the color indicates the magnitude of the attenuation coefficient).

Figure 10 .
Figure 10.Comparison between negative capacitance circuits and traditional LR circuits.

Figure 10 .
Figure 10.Comparison between negative capacitance circuits and traditional LR circuits.
13 and No. 15 exhibit the same current peak.

Table 1 .
Parameters of the substrate plate.